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We follow our general model in Ref. Q and analyze the formation of retinotopic projections for 
the biologically relevant situation of spherical geometries. To this end we elaborate both a linear 
and a nonlinear synergetic analysis which results in order parameter equations for the dynamics of 
connection weights between two spherical cell sheets. We show that these equations of evolution pro- 
vide stable stationary solutions which correspond to retinotopic modes. A further analysis of higher 
modes furnishes proof that our model describes the emergence of a perfect one-to-one retinotopy 
between two spheres. 

PACS numbers: 05.45.-a, 87.18. Hf, 89.75.Fb 



I. INTRODUCTION 

An essential precondition for a correct operation of the nervous system consists in well-ordered neural connections 
between different cell sheets. An example, which has been explored both experimentally and theoretically in detail, 
is the formation of ordered projections between retina and tectum, a part of the brain which plays an important role 
in processing optical information At an initial stage of ontogenesis, retinal ganglion cells have random synaptic 
contacts with the tectum. In the adult animal, however, a so-called retinotopic projection is realized: Neighboring 
cells of the retina project onto neighboring cells of the tectum. A detailed analytical treatment of Haussler and von 
der Malsburg described these ontogenetic processes in terms of self-organization Q . In that work retina and tectum 
were treated as one-dimensional discrete cell arrays. The dynamics of the connection weights between retina and 
tectum were assumed to be governed by the so-called Haussler equations. In Ref. Q we generalized these equations 
of evolution to continuous manifolds of arbitrary geometry and dimension. Furthermore, we performed an extensive 
synergetic analysis 0, Q near the instability of stationary uniform connection weights between retina and tectum. 
The resulting generic order parameter equations served as a starting point for analyzing retinotopic projections 
between Euclidean manifolds in Ref. |i6|. Our results for strings turned out to be analogous to those for discrete 
linear chains, i.e. our model included the special case of Haussler and von der Malsburg Additionally, we could 
show in the case of planar geometries that superimposing two modes under suitable conditions provides a state with 
a pronounced retinotopic character. 

In this paper we apply our general model |^ again to projections between two-dimensional manifolds. Now, however, 
we consider manifolds with constant positive curvature. Typically, the retina represents approximately a hemisphere, 
whereas the tectum has an oval form jj]. Thus, it is biologically reasonable to model both cell sheets by spherical 
manifolds. Without loss of generality we assume that the two cell sheets for retina and tectum are represented by the 
• ^ surfaces of two unit spheres, respectively. Thus, in our model, the corresponding continuously distributed cells are 
represented by unit vectors f and i. Every ordered pair (f, f) is connected by a positive connection weight w{i,f) as 
H , is illustrated in Figure ^ The generalized Haussler equations of Ref. 0, for these connection weights are specified 
as follows 

w{i,f)^f{i,f,w)-^^^ [dnt,f{i',f,w)-'^^ [d^ir'f(t,f',w). (1) 
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The first term on the right-hand side describes cooperative synaptic growth processes, and the other terms stand for 
corresponding competitive growth processes. The total growth rate is defined by 

f{i,r,w) =a + wii,f) Jdnt' Jdnr'CTii-i')cRir-r')wii',f'), (2) 

where a denotes the global growth rate of new synapses onto the tectum, and is the control parameter of our system. 
The cooperativity functions exit ■ t'), c_R(f • r') represent the neural connectivity within each manifold. They are 
assumed to be positive, symmetric with respect to their arguments, and normalized. The integrations in ^ and 
(0) are performed over all points t,f on the manifolds, where d^lt,d^lr represent the differential solid angles of the 
corresponding imit spheres. Note that the factors Stt in Eq. are twice the measure M of the unit sphere, which is 
given by 

2W TT 

Af = jdflt^ jdflr^ jdip JsmM^^AiT. (3) 



If the global growth rate of new synapses onto the tectum a is large enough, the long-time dynamics is determined by 
a uniform connection weight. However, we shall see within a linear analysis in Section ITll that this stationary solution 
becomes unstable at a critical value of the global growth rate. Therefore, we have to perform a nonlinear synergetic 
analysis, in Section IlIII which yields the underlying order parameter equations in the vicinity of this bifurcation. As 
in the case of Euclidean manifolds, we show that they have no quadratic terms, represent a potential dynamics, and 
allow for retinotopic modes. In Section llVI we include the influence of higher modes upon the connection weights, 
which leads to recursion relations for the corresponding amplitudes. If we restrict ourselves to special cooperativity 
functions, the resulting recursion relations can be solved analytically by using the method of generating functions. As 
a result of our analysis we obtain a perfect one-to-one retinotopy if the global growth rate a is decreased to zero. 



II. LINEAR ANALYSIS 

According to the general reasoning in Ref. Q we start with fixing the metric on the manifolds and determine the 
eigenf unctions of the corresponding Laplace-Beltrami operator. Afterwards, we expand the cooperativity functions 
with respect to these eigenf unctions and perform a linear analysis of the stationary uniform state. 



A. Laplace-Beltrami Operator 



For the time being we neglect the distinction between retina and tectum, because the following considera- 
tions are valid for both manifolds. Using spherical coordinates, we write the unit vector on the sphere as 
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X = (sin-iJcost/j, sin?9sin((9, COS"!?). The Laplace-Beltrami operator on a manifold reads quite generally [g 

A = ^9,(/'^V^a^) . (4) 

For the sphere the components of the covariant tensor g/^i, are 

f dx\^ ^ dx dx ^ f dx\'^ . ^ 

311= (^^j =1. .gi2=52i = ^^=0, 522= (^^j =sm^. (5) 

With this the determinant of the covariant metric tensor reads g — sir? d and the components of the contravariant 
metric are given by 

1 

whence the Laplace-Beltrami operator for the sphere takes the well-known form 

sh[d dd y"^"^ ddj^ sin2 {} 

Its eigenfunctions are known to be given by spherical harmonics Y/m (i) : 

^i,,^Yi^{x) = -l{l + \)Yi^{x). (8) 

With Z = 0, 1, 2, . . . and m = — — / + 1, . . . , ^ — 1, ^ they are {21 -f l)-fold degenerate and form a complete orthonormal 
system on the unit sphere: 



11 1 12 21 n 22 ^ (a\ 



ldnxYhn{x)Yi!'^,{x) = Sll'Srnm' , (9) 
oo / 

J2 E yim{£)Yr^{£') - six-x'). m 

1=0 m=-l 

B. Cooperativity Functions 

The argument of the cooperativity functions c{x ■ x') is the scalar product x ■ x' which takes values between —1 and 
+1. Therefore the cooperativity functions can be expanded in terms of Legendre functions Pi{x ■ x'), which form a 
complete orthogonal system on this interval P, 7.221.1]: 

1 

Pi{a)P,{a)da = , (11) 

-1 



^ OO 

-Y,{2l + m{a')Pi{a) = S{a-a'). (12) 

(=0 

Then the expansion of the cooperativity functions read 



°° o; I 1 

c(£-£')=E^^/'^'(*-^')' (13) 



47r 

1=0 



where fi denote the respective expansion coefficients. Using the Legendre addition theorem [T^ 

E Y,^ix)Yr^ix'), (14) 



rn— — l 

we arrive, for each manifold, at the expansion 

CO L oo I 

L=OM=-L 1=0 m=- 1 

Note that the normalization of the cooperativity functions and the orthonormality relations lO lead to the constraints 
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C. Eigenvalues 



The initial state of ontogenesis with randomly distributed synaptic contacts is described by the stationary uniform 
solution of the generalized Haussler equations, WQ{t,r) = 1. Its stability is analyzed by linearizing the Haussler 
equations with respect to the deviation u(i,f) — w{t,r) — WQ{t,f). The resulting linearized equations read 



with the linear operator 

L(t, f, v) = —av{t, f) 



dnt, 



v{t, f) — L{t, f, v) 



v{i',f)+ dQt" dVlr.. CT{i' ■i")cR(f ■ f")v(i" ,f") 



(16) 



/ dVlr' 



v{i, f') + / dnt" / dnr" CT{i ■ i") CR{f' ■ f") v{i",f") 



+ / dVLf / dVLr' CT{i ■ i') CR{f ■ f) v{i', f') . (17) 



To solve Eq. H16|l . we have to consider the eigenvalue problem of the linear operator IjlTfl . It has the eigenfunctions 

vt'r{i,r)^Yl^S)yi^a{r) (18) 

and the spectrum of eigenvalues reads 0: 



A Mm 
LI 



-a -I L = M ^l = m = 

+ iflfP - l)/2 L = M = 0, (l, m) + (0, 0) 
; = m = 0, (L, M) ^ (0, 0) 
— a + fj^fl^ otherwise . 



(19) 



By changing the uniform growth rate a in a suitable way, the real parts of some eigenvalues 1)19(1 become positive and 
the system can be driven to the neighborhood of an instability. Which eigenvalues H19|) become unstable in general 
depends on the respective values of the given expansion coefficients /J, /;^. If we assume monotonically decreasing 
expansion coefficients fj^, 



l = .fI>fI>fI>--->0, 1 = > /f > > • • • > , 



(20) 



the maximum eigenvalue in (|19|l is given by Amax = A]*(™ — —a + fff^ . Thus, the instability occurs when the global 
growth rate reaches its critical value ac = At this instability point all nine modes with = (1, 1) and 

M" = 0, ±1, = 0, ±1 become unstable, where we have introduced the index u for the unstable modes. 



III. NONLINEAR ANALYSIS 

In this section we specialize the generic order parameter equations of Ref. to unit spheres. We observe that 
the quadratic term vanishes and derive selection rules for the appearance of cubic terms. Furthermore, we essentially 
simplify the calculation of the order parameter equations by taking into account the symmetry properties of the cubic 
terms. We show that the order parameter equations represent a potential dynamics, and determine the underlying 
potential. 



A. General Structure of Order Parameter Equations 

The linear stability analysis motivates treating the nonlinear Haussler equations near the instability by decom- 
posing the deviation v{t,r) = w{t^f) — WQ{t,r) in unstable and stable contributions, 

v{if) = U{if) + S{if). (21) 

Using Einstein's sum convention the expansion of the unstable modes reads 

C/(t,f) = U,'r^''Y,l,4i)Y,l4f), (22) 

and, correspondingly, the contribution of the stable modes is given by 

Sii,r)^Sll"-Ylj,,ii)Y,^ir). (23) 
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Note that the summation in H23() is performed over all parameters {L,l) except for (L",?") = (1, 1), i.e. from now 
on the parameters {L,l) stand for the stable modes alone. With the help of the slaving principle of synergetics 
0, 13 the original high-dimensional system can be reduced to a low-dimensional one which only contains the unstable 
amplitudes. The resulting order parameter equations read Q 



U 



AU 



Am" ,m"'m^" jjM^'m"' jjM"^" m"" , rim" ,m"' m"" m"'" rrM"'m"' jjM""m"" jtM'-"" m"'" 
^A/f" /\/f"'A/f"" ^ "T -Ojvf" M"'M""Af ""' ^ ^ 



(24) 



They contain, as usual, a linear, a quadratic, and a cubic term of the order parameters. The corresponding coefficients 



can be expressed in terms of the expansion coefficients fj^, of the cooperativity functions H15|l and integrals over 



products of the eigenfunctions Yim{x): 



•^/(i)i<2)...;( 

The quadratic coefficients read 



d^x Yi*„-^{x) Yi(i) „(i) (x) Yi(2) „(2) (£)••• Yl(,^) „(„> (i) , 
d^x Yi(i) ,„(i) (x) Yii2) „j(2) (x) • • • Yi(^) ,„(„) (x) . 



(25) 
(26) 



A 



m ,m m 



T /-i? tM" ,M'^' M"" Tr?i",m"'m"" 



'M"" ~ Jl Jl -'1,1 1 



1,1 1 



whereas the cubic coefficients are 



(27) 



T^m ,m m m 
i'M",M«"M""M""' 



1 

8^ 



* fT fR f TAr.M"'M""AP"" r Tm""m" 
■^7l Jl 1^1,111 Omi'm'" 



m .ni m m 



+ ^1.1 1 1 



M""J\/""' 



fT i-_R , rT rBl tM" ,M"' M j-m'',m"'m 1 



\ fT fR , /-T /-i?, 1 tM"-,ivi ivi jn 
[JL Jl + /l /l J -'l,lL -'1,1; 



(1 + /I) /^lY^"'^^ 



77 



40F 



R\ jm",rn"'m 
I 



LI 



(28) 



Note that Eq. (|28() involves a summation over all stable modes {L,M;l,m). As is common in synergetics, the cubic 
coefficients (|28|l consist in general of two parts, one stemming from the order parameters themselves and the other 
representing the influence of the center manifold H on the order parameter dynamics according to 



oMm _ jjMm,M"m" M"'m"' tjM" m" jjM" ' r 
LI LI 



(29) 



Here the center manifold coefficients m m m defined by 



H 



LI 



'L,l 1 



2A - An 



'/,1 1 



"^1 



LI 1 



'L,l 1 



(30) 



B. Integrals 

The order parameter equations contain the following integrals: JJ^ ™ , ™ , ™ , 7™j|™ ™ , /™xii ™ ™ • 
The first integral is obtained by the orthonormality relation Q and 

y,_™(i) = (-i)'"r,;,(x), (ai) 

yielding J"l ™ = (— l)™^m',-m" ■ Integrals over three and four spherical harmonics can be calculated with the help 
of the following relation : 



h+h h 



Yi,^^AS')Yi,,^,{x) = J2 E \I^^^^^^^^^C{lu0j2,0\h,0)C{h,mij2,m2\h,m3)Y,^,^,{x), (32) 
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where C(Zi, mi, ^2, ma) represent the Clebsch-Gordan coefBcients IqI. Applying (|32ll to integrals over three 

spherical harmonics leads to 



For V = 1" = 1 it follows 

/™f '™" = ^^^^^^ ^ 0, 1, 0|/, 0) m', 1, m"|Z, m) . (34) 

As the Clebsch-Gordan coefficients C(/i, 0, ^2, OI/3, 0) vanish if the sum I1 + I2 + 13 is odd 0|, we obtain ™ = 0. 
Thus, the quadratic contribution H27|) to the order parameter equations l|24|l vanishes, by analogy with Euclidean 
manifolds 6j. Furthermore, non- vanishing integrals (|34|) can only occur for I — and / = 2. For I = we obtain from 
the Clebsch-Gordan coefficients J^J the result 

jO,m m V ^) sr f'iK\ 

-'0,11 ~ 0„i',m" ■ ['J'i) 



For I = 2 we find, correspondingly, the nonvanishing integrals 



rO,00 _ 1 tO,00 _ 1 tO,1-1 _ rO,-ll _ 1 rO,l-l _ rO,-ll 

^0,11 ~ /I— ' -'2.11 ~ rr- ' 0.11 ~ ^0,11 ~ /I— ' ^2,11 ~ -'2,11 



rl,10 _ rl,01 _ T-1,-10 _ T-1,0-1 

-'2,11 — ^241 — ''2,11 — ''2,11 



3 t2,11 _ 3 r-2.-l-l 

7 ^2.11 — 



2V157r 



/SOtt 



1 -'2,11 



/SOtt 



(36) 



Furthermore, the integrals Ij^ follow from 



Integrals over four spherical harmonics can also be calculated with the help of H32|l . and the result is 

(3 



E E 

is — ll" —I'" \ rn3 — — l3 



^ (21" + i){2r" + 1) 

4^(2^3 + 1) 



C{1", 0, 0|?3, 0)C{1", m", Z'", m"'|Z3, ma)/™'^ 



(37) 



(38) 



Specialyzing to Z = /' = /" = /"' = 1 and taking into account leads to /{"i™/" " oc (5m'+m"+m"',m- Thus, 

we obtain the selection rule that the nonvanishing integrals fulfill the condition m' + m" + m'" — m. The 

detailed evaluation yields for those the respective values 



rO.OOO 

^1,111 

tO,1~10 _ tO,-110 _ tO,10-1 _ xO,-101 _ xO,01-l _ rO,0-ll 
-"l,!!! ~ -'1,111 ~ -"1,111 ~ -'1,111 ~ -"1,111 ^ -'1,111 

rl.lOO _ t1,010 _ tI.OOI _ T~l,-100 _ T-1,0-10 _ r-1,00-1 

-"1,111 — -"1,111 — -"1,111 — -'1,111 ^ -'1,111 — -"1,111 

t1,11-1 _ rl,l-ll _ rl,-lll _ x-1, 1-1-1 _ T-1,-11-1 _ T-1,-1-11 

-"1,111 — -"1,111 — -"1,111 — -"1.111 — -"1,111 — -"1.111 
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20^ ' 
3 

^20^' 
3 

20^ ' 
3 



(39) 



C. Order Parameter Equations 

To simplify the calculation of the cubic coefficients 1)28(1 in the order parameter equations ((24|l . we perform some 
basic considerations which lead to helpful symmetry properties. To this end we start with replacing by — m". Using 

It u' u" u'" u 11' u" u'" 

Eq. (|31|l we obtain /™-^ ^™ " ™ — IiTi i~™ ~™ • Corresponding symmetry relations can also be derived 

for the other terms in H28|) . Therefore, we conclude that the order parameter equation for {/-^f"-™" jg obtained from 
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that of by negating all indices M" and m" with unchanged factors. Thus, instead of explicitly calculating 

nine order parameter equations, it is sufficient to restrict oneself determining the order parameter equations for U'^^, 
Ijw^ IjOi^ Q^j^^ ljii_ Xhe remaining five order parameter equations follow instantaneously from those by applying the 
symmetry relations. With this the order parameter equations result in 

jjoo ^ ^^00 ^ ^^([/00)3 _ 2/32C/™C/-i°C/^° - 2/32 C/"° + 2P3U°°U^-'U~^' + 2/33C/™C/-^-^C/" 
[>" = AC/"+/34C/°°C/°^C/^" + /35(C/"')'C/'-'+/?6C/"iC/"^'C/"+/93(C/°°fC/" 

+/35(f/-'°)'C/'-' + ^6(7-i";7i";7-i-i + PrU-'-'U-^'U'-' + ps{U-'-'fU'' , 
[ji-i = Af/i-i + /34C/°"C/"-iC/i" + p5{U°-')^U'' + PeU°-'U°'U'-' + P^iU^'^fU^-^ 
+/35(C/i°)'C/-i-i + /36C^'"t^"'"f/'"' + PiU^-^U^^U-^-^ + /?8(C/i-i)'C/-" , 
= A[/-" + /34C/*'°C/*^iC/-i° + p5{U°^fU-^-^ + /^gL/Oif/O-iL/-" + /^al;/"")^;/-" 
+/35(C/-^")'C/" + /36(7-i"(7i";7-ii + /37[/-"[/-i-i[/" + P^iU-^^U^^^ , 

f7"i = A;7°i + ^2C/"'(c/"°)' + /39(c/"')'f/""'-2/33;7°i;7i";7-i"-/34;7°"[/"[/-i° 

-/34{/"°{/i°{/-" - peU"'U''U''-' - PgU'^'U'-'U-'' - 2(3,U"-'U''U-'' , 
^0-1 ^ A[/"-i + ;32C/""'(C/"°)' + /39(C/*'-i)2C/*'i - 2/33;7°-i;7-i"f/i" - /34C/°°f/-i-iC/i" 

-piU°"u-'"u'-' - /36;7°-^f/-^-^f/" - /36(7"-i;7-";7i-i - 2PzU''^u-^-^u^-^ , 

f/10 = A[/i"+/32C/i°(C/*'°)'+;59(C/'")'C/"'""2/33t/i";7"i;7"-i-/34f/°*'f/"f/"-' 

-/34C/°°C/°^C/^-^ - PeU'^'U^'U-'-' - /36(7i"t/-"!7i-i - 2/35C/-^°C/"C/^-^ , 
f7-i" = A[/-i" + /32C/-i°(C/*'°)' + ;99([/-i°)2[/io - 2/33[/-i°f/"-i[/"i ~ /34C/*'"[/-i-iC/°i 

-/34f/°°t/°-'C/-" - ^6(7-i"[/-i-iC/" - ^eLZ-i^L/i-iL/-" - 2p,U''U-'-'U-'' . (40) 
With the abbreviations 7 = j/'^''^, 7 — ff fi ^^'^ 7^'' — flfP respective coefficients in (|40ll read 
9. 7 27 + 72^0-1 7 27 + 70'2_i ^ -.^^2,2 

^ ~80^ ' 

= "80^ 
9 

P^ = -^7 + 



80 27 


— a 


_(^2,0_l)/2 + 


80 27 


-a-(70.2-l)/2 


25 27- 


2 2 

a — 'Y^'^ 


7 


27 


+ 72^°-l 


7 


27 + 70^2 _ 1 


27 7 


+ ^2,2 


40 27 


— a 


- (7^'" - l)/2 + 


80 27 


-a-(70.2_i)/2 


25 27- 


a — 72'2 


7 


27 


+ 72^0-1 


7 


27 + 70'2 _ 1 


27 7 


+ ^2,2 



= -^7 



80 ' ' 8027-a- (72^0- l)/2 40 27-a-(70'2_l)/2 25 27 - a - 72.2 ' 
Q ^ 27 + 72^0 -1 7 27 + 7°'2_i 11-^ 7 + 72-2 



80' 16027-a-(72.o_i)/2 160 27 - a - (7O.2 _ i)/2 20027-0-7 

3. 37 27 + 72.0 -1 37 27 + 70^2 _i 2I7 7 + 72-2 

Pi = ~77T7 



2,2 



40' I6O27- a -(72.0-1)72 160 27-a- (70.2_i)/2 20027-0-72.2 

3 , 37 27 + 72-0 - 1 37 27 + 7".2 _ 1 3^ ^ ^ ^2,2 

P5 = 77:7 ' 



40 ' 160 27 - a - (72.0 - 1) /2 160 27 - a - (70.2 - l)/2 200 27 - a - 72.2 

_ _3_, 7 27 + 72.0-1 7 27 + 7°.2_i i3;y 7 + 72.2 

~ 20^^ 32 27- a -(72.0-1)72 160 27- a- (70.2-1)72 20027-0-72.2' 

5 3,7 27 + 72.0 -1 7 27 + 7O.2-I 137 7 + 72.2 

Pe = 7^7:7 



20' 160 27- a- (72.0 - 1)72 3227- a- (7O.2 _ 1)72 20027-0-72.2 
3. 7 27 + 72.0 -1 7 27 + 7O.2-I 1I7 7 + 72.2 

Pt = -T7T7 



P& = ^7;7t7 



10' 32 27- a- (72.0 - 1)72 32 27- a -(70.2-1)72 20027-0-72.2 

Q 7 27 + 72.0-1 7 27 + 70.2-1 197 7 + 72.2 



20 ' 160 27 - o - (72.0 - 1) 72 160 27 - o - (7O.2 - 1) 72 200 27 - o - 72.2 
9, 7 27 + 72.0 -1 7 27 + 7O.2-I 7 7 + 72.2 

P9 = 77:7 



40' 40 27- o- (72.0 -1)72 80 27-0 -(70.2-1)72 2527-0-72.2' 
^ 9 7 27 + 72.0-I 7 27 + 7O.2-I 7 7 + 72.2 

40^ 80 27-0 -(72.0 -1)72 40 27-0 -(70.2-1)72 ^25 27- 0-72.2- ^ > 
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The first term proportional to 7 describes the influence of the order parameters themselves, while the other terms 
stand for the contributions of the center manifold. 



D. Real Variables 



To investigate how the complex order parameter equations contribute to the one-to-one retinotopy, we transform 
them to real variables according to 



- U 

- u 



uo = C/"7V2 , Ul - ((7" + C/-i-i)/2 , U2 

UQ = i{U"^ + U°-^)/2 , U7 = (L/i"-[/-i")/2 , Us 
Then the equations of evolution for the real variables Ui read 

Uo = Auq + 2f3iul + 2(32Uo{ul + uj) + 2/32Mo(w? + ul) + 2f33Uo{ul + ul + u^ + u\ 

+ \/2 /34(miU5U7 + + U2UqUj + U^UqUj — UiUqUs — U^U^Ur — U4U5U8 - 

Ul = Aui + V2 Piuoiusur - uqUs) + P^iusul - Usui - 2u4U5Ue) - Peuiiul + ul) 
+2P3ului + Pziu^u^ - Usui + 2M4M7U8) + (37Ui{ul + ul) + /3sui{ul + uj) , 

U2 = Au2 + V2 PiUo{u5Us + uqUy) + (i^^ului - Uiul + 2u3U5Ue) - l3eU2{ul + ul) 
+2P3ulu2 - Pziu^u^ - Uiul - 2u3U7Us) + l37U2iul + ul) + f3sU2{ul + ul) , 

U3 = Au3 - V2 PiUoiusUr + uqUs) + I35{uiul - uiul + 2u2U5Ue) - Peusiul + ul) 

+ 2/33MqU3 -I- Psiuiu"^ - Ulul + 2u2U7Us) + (i7U3{u\ + ul) + Ps,U3{ul + ul) , 

U4 = Aui + V2 PiUo{uQUT - u^us) + I3z{u2ul - U2ul - 2uiUr^UQ) - l3iiU4,{ul + u^) 
-|-2/33MqU4 -I- (iz{u2ul — U2M7 + 2uiU7Us) + P7U4{ul + ul) + PsU4{ul + ul) , 

U5 = Au5 + 2/32Wo"5 - Pgusiul + ul) + 2/33U5(uy + ul) + V2 /3iUQ{uiU7 - usur + U2Us 
-PqU5{uI + ul + ul + ul) - 2[3z{uiUaUq - U1M3U5 - U2U4W5 - , 



U6 = Aug + 2/32MoU6 - /39U6(U5 + We) + 2/33U6(w7 + Ws) + ^/2/34Uo(u2W7 ■ 



U4U7 ■ 



~PqUq{ui -f U2 + "3 + "4) ^ 2(3^{uiUaUz + U1U3UQ + U2UAUQ - U2U3U5) , 
iij — Au7 + 2l32uluj — PqUt^u^ + ul) + 2l33Uj{ul + ul) + V2 l34Ua{uiU5 — U3U 

—(3eU7{ul + ul + ul + ul) — 2(}^{u2UiU7 — U1U3U7 — U2U3US ~ uxu^u^) , 
Us = Au8 + 2/32UoUs-/39Us(u7 + u|) + 2/33Us(u^+u^)- V2/?4Uo(uiU6-l-U3U 

-PqUs{uI -f U2 -f U3 -I- ul) - 2/35(uiU3Us - U2M3U7 - U1U4M7 - U2U4U8) . 

Note that the real order parameter equations 1021 ISJ follow according to 

. _ dV{{u,}) 



from the potential 
Vi{u,}) ^ 



A 



/3i 



dui 



4 5 2/2, 2\ a 2/ 2 1 2\ a 2/ 2 , 2 , 2 , 2\ 
io - P2Uq[u^ + Ug) - P2Uq{Uj + Ug) - /^3Uo("l + "2 + ""3 + W4) 



-i-i)/2 




-i)/2 


(42) 


-^")/2 




:) 

- U3U6M8) , 


(43) 


- /36Ui(uy -t- ul) 






(44) 


- /36U2(Uy -f u|) 






(45) 


- /36U3(uy + u|) 






(46) 


- [3(iUi{u^ -f Ug) 






(47) 


- U2U8 - U4U8) 






(48) 


- uius - U3U8) 






(49) 


- U2U6 + U4U6j 






(50) 


- U2U5 -f U4U5) 






(51) 




(52) 






\^ul) 





-\/2j3iUQ{uiUzU7 + U2U5M8 + U2UQU7 + U4V.QU7 — UiUgUs — U3U5U7 — U4U5U8 — U3UgM8) 

-I3z{ul - Ug)(uiU3 -I- U2U4) - (3z{uj - ul){uiU3 - U2U4) - 2/35U7U8(uiU4 + U2U3) 



-2/35U5Ue(u2U3 - U1U4) -t- -[/3g(u5 -I- Ug) -f /3e(u7 + U8)](ui -f- U2 



Pi, 2 



ul){u 



o / 2 I 2\i 2 I 2\ \ / 2 , 2\2 , / 2 , 2\21 , / 2 , 2\2 , / 2 , 2n2 

-/33(% +U6)(U7 +U8) - — [(Ui +U2) + (Ug +U4) J + — (Ug +Ug) +— (Uy + Us) 



(53) 



Naturally, a complete analytical determination of all stationary states of the real order parameter equations H43|l - H51() 
is impossible. However, we are able to demonstrate that certain stationary states admit for retinotopic modes. 
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E. Special Case 

To this end we consider the special case ui,U2,U5,uq,U7,us — 0. Then the equations (03, and l|T7|l for the 

non- vanishing amphtudes uq, "3, U4 reduce to 

uo = Auo + 2f3iul + 2f33{ul+ul)uo, 
U3 = Au3 + 2P3ulu3 + f3siul + ufjus , 

U4 = Au4 + 2f33ulu4 + Ps{ul + 'U'l)u4 . (54) 

Due to the relation 

^ = ^ (55) 

U3 U4 

one obtains constant phase-shift angles, i.e. it holds U3 oc U4. Therefore, the system of three coupled differential 
equations can be reduced to two variables. To this end we introduce the new variable 



^^^uj + ul, (56) 

which leads to 

Uo = Auo + 2f3iul + 2M'^uo , 

e = AC + 2/33ulC + M'- (57) 

The stationary solution, which corresponds to a coexistence of the two modes, is given by 

2 _ A .2 _ A 

2(/33 + /?8)' ^ " /33 + /38' ^^^> 

where we used the relation /Jg = Pi +03 following from (|41ll . Demanding real amplitudes uq, ^ leads to the coexistence 
condition 

/33 + /38 < . (59) 

Furthermore, we require stability for this state. Therefore we consider the corresponding potential V{uq,^), which 
can be read off from (|55)l and l|5B|) : 

v{uo, o--^{ui+ e) - ~ Psuie ~ . (eo) 

Stable states correspond to a minimum of V, which leads to the conditions 

2/33 - /38 > , /33 - A > . (61) 

The inequalities H59() . H61|l can be summarized according to 

/38<0, P3<-P8, 2/33 >/38. (62) 

If they are valid, both the uq- and the f-mode coexist. If we set U4 = 0, without loss of generality, the solution reads 
in complex variables according to (|55|l 



Thus, the unstable part is given by 



uiif) = \yio{i)y^,{f) ~ yii{i)yi-i{f) - yi-i{i)yii{f) ■ m 



P3 + Ps 

Using the Legendre addition theorem H14|l reduces l|64(l to 



U{if)^J--^Piii-f) (65) 
y P3 + P8 

with Pi (t-r) = t-f. Thus, the unstable part is minimal, if t and f are antiparallel, i.e. the distance of the corresponding 
points on the unit sphere is maximum. Decreasing of the angle between t and f leads to increasing values of U{t, f), 
and the maximum occurs for parallel unit vectors. This justifies calling the mode (|65f) retinotopic. 
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IV. ONE-TO-ONE RETINOTOPY 

Now we investigate whether the generahzed Haussler equations 1^ describe the emergence of a perfect one-to-one 
retinotopy between two spheres. To this end we follow the unpublished suggestions of Ref. 0| and treat systematically 
the contribution of higher modes. Because the Legendre functions form a complete orthogonal system (|ll|l . (|12|l for 
functions defined on the interval [— 1,+1], their products can always be written as linear combinations of Legendre 
functions. This motivates that the influence of higher modes upon the connection weights, which obey the generalized 
Haussler equations can be included by the ansatz 



(66) 



1=0 



where the amplitudes Zi are time dependent. 



A. Recursion Relations 

Inserting (|66|) into the generalized Haussler equations and performing the integrals over the respective unit 
spheres leads to 



^(2? + l)Z,P,(a) =a 



1=0 



l-J2{2l + l)ZiPi{a) 



Y^i^l + \)ZiPi{a) + l)Z,f^f^\P,{a) - Z,] . (67) 



1=0 



l'=0 



The products of Legendre functions occuring in (|67|l can be reduced to linear combinations of single Legendre hmctions 
according to the standard decomposition 9, 8.915] 



Pi{a)Pv{a) =J2^iM,kPi+i'-2ki<y), I < I' 



k=0 



with the coefficients 



{21' + 2l-4k + l)ai,^kakai^k (2fc-l)!! 

^1,1', k — /„„ , — T , Ctfe — 



(2;' + 2/-2fc+l)ai+i/ 



fc! 



(68) 



(69) 



Thus, contributions to the polynomial Pi^o) only occur iff the relation fc = {I + V — l)/2 is fulfilled. Furthermore, 
using the orthonormality relation Hll|l yields the following recursion relation for the amplitudes Zi : 



{21 + l)Zi = a[Si^„ - {21 + l)Zi\ - {21 + l)Zi{Zl + Sf^f^Zf) + ^(2/' + l)Zi, 



l'=0 



^(2r + i)z,,/,^/« 



U"=0 



k=0 



l"=l'+l 



k=0 



(70) 



Note that Eq. H70(l cannot be solved analytically for arbitrary expansion coefRcients ff, ff^ of the cooperativity 
functions. Therefore, we restrict ourselves from now on to a special case. 



B. Special Cooperativity Functions 

For simplicity we assume that the expansion of the cooperativity functions ()13fl breaks down after the first order: 

(71) 

4:71 

With this choice the recursion relation H70|) for I = reduces to 

Zo^-{a + Z'^ + 3-/Z^){Za-l), (72) 



Mi ■ i') = ^[1 + 3/rPi(t • i')] , CR{f .f') = ^[l + Zf?P^{f ■ f')] 
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where we have used again the abbreviation 7 = ff fi . For / 7^ 0, by taking into account H69|) . we obtain 

^ -{a + Zl + i^ZDZi + Z,Zi + 37Z1 ^^'''^2/+"^/^^'^' ■ (^^^ 

The long-time behavior of the system corresponds to its stationary states. They are determined hy Zq = 1 from H72|l . 
whereas ()73|l leads to a nonhnear recursion relation for the amplitudes Zi with I =^ Q. However, by introducing the 
variable 

+ (74) 

this nonlinear recursion relation can be formally transformed into the linear one 

{l + l)Zi+^{u) = {2l + l)uZi{u) + lZi^i{u), 1>1. (75) 

Thus, solving the nonlinear recursion relation H73|l amounts to solving the linear recursion relation H75f) for Zi{u) in 
such a way that the self-consistency condition H74(l is fulfilled. 



C. Generating Function 

To determine the amplitudes Zi (u) we calculate their generating function 

E{x,u) = "^Zi {u)x' , (76) 

where we have the normalization 

E{0,u) = Zoiu) = l. (77) 

Multiplying both sides of H75|l with and summing over / > 1 leads to an inhomogeneous nonlinear partial differential 
equation of first order for the generating function: 

(x2 - 2ux + 1) ^^1^ = _ x)E{x, u) + Zi{u)-u. (78) 
ox 

At first, we consider the homogeneous equation corresponding to H78I) : 

{x' ~ 2ux + 1) ^ - x)E^,^{x, u) . (79) 

It is solved by the method of separating variables, yielding 

K(u) 

Ei,,„,{x, u) - \' , (80) 

yx — Zux + 1 

where K{u) is an integration constant. Afterwards, we determine a particular solution of the inhomogeneous equation 
(|78ll by using the method of varying constants. Using the ansatz 

„ , N K(x,u) 

E^,,,{x,u)^^ \ ' (81) 
va:^ — 2ux + 1 



leads to the differential equation 



dKjx^u) ^ Zi{u) - u 
dx y/x^ -2ux + l' 



which is solved by using !?, 2.261]: 



[Ziju) - u] \n[2Vx^ -2UX + 1 + 2{x - u)] 
""("'"^^ Vx^-2.x+l ■ ^''^ 
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Thus, the complete solution E{x, u) = Eho^{x, u) + £'pait(a;, u) of Eq. H78() reads as follows: 



^, , K{u) + [Zi{u) - u] ln[2Va;2 - 2ua; + 1 + 2{x - u)] , ^, 

^(2;, u) = ===== . (84) 

Va;^ — 2ux + 1 

Furthermore, using the normalization condition (|77|l fixes the integration constant to K{u) = 1 — [Zi [u) — u] ln(2 — 2u) . 
Thus, the generating function is finally given by 



1 + [Zi [uj ~ u\ In 

E{x,u) = ^ — ^ . (85) 

Vx^ — 2ua; + 1 

D. Decomposition 

We now determine the unknown amplitudes Zi{u). From the mathematical literature it is well-known that the 
recursion relation lf75|) holds both for the Legendre functions of first kind Pi{u) and second kind Qi{u), respectively 
Thus, we expect that the generating function H85() can be represented as a linear combination of the generating 
functions of the Legendre functions of both first and second kind, which are given by j^, 8.921] and |9|, 8.791.2]: 

00 ^ 

Ep{x,u) = ^Pi(^)x' = ^ (86) 
^ Va; — 2ux + 1 

\J x^ — 2ux + 1 + u — X 
00 In — 



V a;^ — 2ux + 1 



Indeed, taking into account the explicit form of the Legendre function of second kind for ^ = [ll 

, ^ 1 u+1 
^°(")=2^"^' 

the generating function (|85|) decomposes according to 

E{x,u) = {l + [Zi{u) - u]Qoiu)}Ep{x,u) - [Zi{u) - u]EQix,u) . (89) 

Inserting H86|) , H87|l and performing a comparison with (|76|l then yields the result 

Zi{u) = {l + [Zi{u) - u]Qoiu)}Pi{u) - [Ziiu) - u]Qiiu) . (90) 

Thus, the amplitudes Zi{u) turn out to be linear combinations of Pi{u) and Qi{u). To fix the yet undetermined 
amplitude Zi[u) in the expansion coefficients of H90() . we have to take into account the boundary condition that the 
sum in the ansatz has to converge. 

E. Boundary Condition 

Because the Legendre functions Pi{<7) do not vanish with increasing I, we must require 

limZ,(u) = 0. (91) 

I — >oo 

The series of Legendre functions of first kind Pi {u) with fixed u > 1 diverges for / 00 according to 0, 8.917] 

Po(u) <Pi(w) <P2(m) < ... <P«(m) < ... , u>l. (92) 

The Legendre functions of second kind Qi{u), however, converge to zero (see Figure |2Il. Thus, performing the limit 
Z ^ 00 in Eq. (|^ . we obtain 

1 + [Zi(u)-m]QoH = 0. (93) 
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FIG. 2: The Legendre functions of first and second kind Pi{u) and Qi(u) for u> 1. We have Pi(l) = 1, whereas Qi{u) diverges 
for u [ 1. Important for the boundary condition of Zi(u) is the different behavior for increasing values of I: Piiu) diverges 
according to H92|l . whereas Qi(u) converges to zero. 



From the explicit form [13j Qi{u) — uQq{u) — 1 it follows that Zi{u) is fixed according to 

Qi{u) 



Zi{u) = 



Qoiu) 



With this we obtain that the result finally reads 



Zi{u) 



Qi{u) 
Qo{u) 



which is not valid only for / 7^ but also for / = due to 177|l . 



(94) 



(95) 



F. Connection Weight 
Inserting (|95|l into (|66|l yields the following solution for the connection weight: 

Qo(u) ^ 



1=0 



Using the identity J„ 8.791.1] 



00 ^ 

yi2l + l)Qi{u)Pi{<j) = 



and H88|l . we obtain for the connection weight 



w{a) 



2 A u + 1 
In ■ 



M — (7 V U — 1 

Note that integrating H98|) over the unit sphere leads to 

27r +1 

J'^'^ j'^^ w(cr) = 47r , 
-1 



(96) 



(97) 



(98) 



(99) 



i.e. the total connection weight coincides with the measure Q. 
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FIG. 3: a) Relation 1)100^ between the control parameter a and the variable u. b) The connection weight for diflFerent values 
of the control parameter a. For decreasing values of a the connection weight around a — +1 is growing. In the limiting case 
a —> the connection weight w{a) becomes Dirac's delta function Hl()7|l . 



On the other hand we have to take into account that the self-consistency condition (|74|l yields an explicit relation 
between the variable u and the control parameter a. Indeed, we infer from (|74|l and 1)94(1 the following transcendental 
relation between a and u 



7 3 V " - 1 



which is depicted in Figure 



2 In- 



1 



(100) 



G. Limiting Cases 

The limiting value of (|1UU|I for m — > oo is determined with the help of the expansion j^, 1.513] 

1 -U CO ^ 

k=l 

and reads 

lim a = 7. (102) 

u — ^ 00 

Thus, we conclude that the case u —> 00 corresponds to the instability point ac — /i"//^, which was obtained from 
the linear stability analysis in Section ^1 Correspondingly, using again H100|) . we observe that the connection weight 
()98|l coincides in the limit u 00 with a uniform distribution: 

lim w{a) = 1 . (103) 

Another biological important special case is m J, 1, where we obtain from (|100|l 

lima = 0. (104) 

nil 

Furthermore, considering the limit m J, 1 in H98|l for ct 7^ u, we obtain 

2 / M + r 



lim In = . (105) 

till u — a \ u — 1 



On the other hand, integrating H98|l for it 4 1 over a yields 



1 1 



lim / ( In ^iii 1 da = 2. (106) 

uii J u - a \ u - 1 ' 
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Therefore, we conclude that the connection weight (|98|l becomes in this hmit Dirac's delta function: 

limw(cr) = 4(5(cr- 1) . (107) 

alO 

Thus, decreasing the control parameter a means that the projection between two spheres becomes sharper and sharper 
(see FigurelHt). A perfect one-to-one retinotopy is achieved for a = when the uniform and undifferentiated formation 
of new synapses onto the tectum is completely terminated. 



V. SUMMARY 



In this series of three papers we have analyzed in detail the self-organized formation of retinotopic projections 
between manifolds of different geometries. Applying our generalized Haussler equations 3] to Euclidean manifolds 
0, and to spheres in the present paper, led to remarkably analogous results. Both for one-dimensional strings and for 
spheres we have furnished proof that our generalized Haussler equations describe, indeed, the emergence of a perfect 
one-to-one retinotopy. Furthermore, we have shown in both cases that the underlying order parameter equations follow 
from a potential dynamics and do not contain quadratic terms. However, in contrast to strings, spherical manifolds 
represent a more adequate description for retina and tectum. Therefore, the present paper represents an essential 
progress in the understanding of the ontogenetic development of neural connections between retina and tectum. 
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